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Abstract 

Two different reasons make the search 
for transients in the nearby Universe (d^ 
200 Mpc) interesting and urgent. First, there 
exists a large gap in the luminosity of the 
brightest novae (—10 mag) and that of sub- 
luminous supernovae (—16 mag). However, 
theory and reasonable speculation point to 
several potential classes of objects in this 
"gap". Such objects are best found in the 
Local Universe. Second, the nascent field 
of Gravitational Wave (GW) astronomy and 
the budding fields of Ultra-high energy cos- 
mic rays, TeV photons, astrophysical neutri- 
nos are likewise limited to the Local Universe 
by physical effects (GZK effect, photon pair 
production) or instrumental sensitivity (neu- 
trino and GW). Unfortunately, the localiza- 
tion of these new telescopes is poor and pre- 
cludes identification of the host galaxy (with 
attendant loss of distance and physical diag- 
nostics). Both goals can be met with wide 
field imaging telescopes acting in concert with 
follow-up telescopes. Astronomers must also 
embark upon completing the census of galax- 
ies in the nearby Universe. 

1 A Historical Summary 

Variable objects have played a major role in 
the history of astronomy. At the beginning of 
the previous century, the pulsating Cepheids 
showed that the Galaxy was much larger than 
had been assumed before. Bubble's discovery 
of Cepheids in the Andromeda nebula showed 
that it was a galaxy similar in size to our own. 

F. Zwicky's 18-inch Palomar Schmidt pro- 
gram was the first systematic study of the 
transient sky. Zwicky laid the foundation of 
the field of supernovae which in turn touched 
upon the origin of cosmic rays, the synthe- 
sis of Iron and higher-(y4, Z) elements and 



the collapse of stellar cores. This success led 
to the 48-inch Palomar Oschin Schmidt tele- 
scope, famous for two Sky Surveys. 

At the end of the previous century, su- 
pernovae came back into main stream. The 
first indication of a new constituent of the 
Universe, dark energy, was deduced from the 
dimming of la supernovae located at cosmo- 
logical distances (Perlmutter et al. 1999, 
Riess et al. 1998). The last decade saw a 
veritable explosion in the field of gamma ray 
bursts - the most relativistic explosions in 
Nature. 

The purpose of this white paper is to ex- 
phcitly draw attention to the importance 
of a systematic study of transients in the 
local Universe (defined here as distance ^ 
200 Mpc). The rationale for this myopic view 
is laid in [|2]and ^ We argue that the path 
to discovery necessarily lies through identi- 
fication and elimination of fore- and back- 
ground transient events (Q. We discuss up- 
coming facilities in (§SD- While detection of 
any cosmic source by LIGO constitutes a ma- 
jor step there is wide-spread agreement of the 
necessity of arc-second localization (which 
means in essence detecting an electromag- 
netic counterpart). We end with a proposed 
plan to identify electromagnetic counterparts 
to LIGO sources (©. 

Here, we limit the discussion to transients 
in the local Universe. We refer the reader to 
the following white papers: Cosmological ex- 
plosions (Wozniak et al.), gravitational wave 
sources (Bloom et al., Phinney et al.) and 
Galactic variable stars (Becker et al.). 

2 Rationale & Motiva- 
tion: Events in the Gap 

A plot of the peak luminosity versus a dura- 
tion that is characteristic (based on physics or 
convention) is a convenient way to summarize 
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Figure 1 1 The phase space of cosmic transients : peak V- 
band luminosity as a function of duration, with color a measure 
of the true color at maximum. Shown are the known explo- 
sive (supernovae) and eruptive (novae, luminous blue variables 
[LBV]) transients. Also shown are new types of transients 
(all found over the last two years): the peculiar transients 
M85OT2006-1, M31-RV and V838Mon, which possibly form 
a new class of "luminous red novae" ; almost every scenario has 
been suggested - core collapse, common envelope event, planet 
plunging into star, a peculiar nova and a peculiar AGB phase; 
the baffling transient with a spectrum of a red-shifted Carbon 
star, SCP 06F6 (Barbary et al 2008; see also Soker et al 2008); 
a possible accretion induced collapse event SN 2005E (Perets et 
al. 2009); the extremely faint, possibly Type Ibc, SN2008ha 
(Valenti et al. 2008); and peculiar eruptive events with ex- 
tremely red progenitors SN 2008S and NGC300-OT (Thomp- 
son et al. 2008, Smith et al. 2008, Bond et al. 2009). [Figure 
adapted from Kulkarni et al. (2007).] 



explosive events. We first focus on novae and 
supernovae of type la (SNIa). As can be seen 
from Figure [1] novae and SN la form a dis- 
tinctly different locus. Brighter supernovae 
take a longer time to evolve (the "Phillips" 
relation; Phillips 1993) whereas the opposite 
is true of novae: the faster the nova decays 
the higher the luminosity (the "Maximum 
Magnitude Rate of Decline", MMRD rela- 
tion; see, for example, Delia Valle and Livio 
1995, Downes and Duerbeck 2000). 

The primary physical parameter in SNIa 
is the amount of Nickel that was synthesized. 
There is almost a factor of 10 variation be- 
tween the brightest ("1991T-like") and the 
dimmest ( "1991bg-like") SNIa. The Phillips 



relation has been quantified with high pre- 
cision and the theory is well understood. In 
contrast, the MMRD does not enjoy the same 
quantity or quality of light curves as those 
of la supernovae. Fortunately, dedicated on- 
going nova searches in M31 and the P60- 
FasTING project (see ^ has vastly increased 
the number of well-sampled light curves. 

The large gap between the brightest novae 
(say Mv ~ —10) and the faintest supernovae 
(My ~ —16), especially on timescale shorter 
than thirty days, is hard to miss. The gap 
beckons astronomers to search the heavens for 
new populations. 

A discussion of potential new classes of 
events in the gap would benefit from a review 
of the basic physics of explosions. An impor- 
tant complication is a potential heat source 
at the center: a hot white dwarf (novae) or 
gradual release of radioactive energy (super- 
novae). 

The primary physical parameters are: the 
mass of the ejecta (Mej), the velocity of the 
ejecta (w^), the radius of the progenitor star 
(-Ro) and the total energy of the explosion 
Two distinct sources of energy con- 
tribute to the explosive energy: the kinetic 
energy of the ejecta, E\ = (l/2)Mejff and 
the energy in the photons (at the time of the 
explosion), Sph- 

Assuming spherical symmetry and homo- 
geneous density, the following equation de- 
scribes the gains and losses suffered by the 
store of heat (E): 

E = e{t)Mei - L{t) - AnR{t)^Pv{t). (1) 

Here, L{t) is the luminosity radiated at the 
surface and e{t) is heating rate (energy per 
unit time) per gram from any source of en- 
ergy (e.g. radioactivity or a long-lived central 
source). P is the total pressure and is given 
by the sum of gas and photon pressure. 

Next, we resort to the so-called "diffusion" 
approximation (see Arnett 1996; Padmanab- 
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han 2000, volume 11, 54.81 



(2) 



where E'ph = aT'^V is the energy in photons 
[V is the volume, (47r/3)-R^], and 



BkM^^/cR 



(3) 



is the timescale for a photon to diffuse from 
the center to the surface. The pre-factor B 
in Equation [3] depends on the geometry and, 
following Padmanabhan, we set B = 0.07. k, 
is the mass opacity. 

We will make one simplifying assumption: 
most of the acceleration of the ejecta takes 
place on the initial hydrodynamic timescale, 
Th = Rq/vs, and subsequently coasts at 
R{t) = Ro + v^t. 

First, let us consider a "pure" explosion i.e. 
no subsequent heating [e(t) = 0]. If photon 
pressure dominates then P = 1/3{E/V) and 
an analytical formula for L{t) can be obtained 
(Padmanabhan, op cit): 



L{t) = Lq exp 



tTh + tV2 



(4) 



here, = B{KMej/cRo) is the initial diffu- 
sion timescale and Lq = Sph/ra- 

From Equation m one can see that the light 
curve is divided into (1) a plateau phase 
which lasts until about r = ^Jr^fh after 
which (2) the luminosity undergoes a (faster 
than) exponential decay. The duration of the 
plateau phase is 



(5) 



and is independent of i?o- The plateau lumi- 
nosity is 

= Wr. = ^^- (6) 

As can be seen from Equation [6] the peak lu- 
minosity is independent of the mass of the 



ejecta but directly proportional to Rq. To 
the extent that there is rough equipartitioijf] 
between the kinetic energy and the energy in 
photons, the luminosity is proportional to the 
square of the final coasting speed, f^. 

Pure explosions satisfactorily account for 
supernovae of type Hp. Note that since Lp oc 
Rq the larger the star the higher the peak 
luminosity. SN 2006gy, one of the brightest 
supernovae, can be explained by invoking an 
explosion in a "star" which is much larger 
(160 AU) than any star (likely the material 
shed by a massive star prior to its death; see 
Smith & McCray 2007). 

Conversely, pure explosions resulting from 
the deaths of compact stars (e.g. neutron 
stars, white dwarfs or even stars with ra- 
dius similar to that of the Sun) will be very 
faint. For such progenitors visibility in the 
sky would require some sort of additional sub- 
sequent heat input and discussed next. 

First we will consider "supernova" -like 
events i.e. events in which the resulting de- 
bris is heated by radioactivity. One can easily 
imagine a continuation of the la supernova se- 
quence. We consider three possible examples 
for which we expect a smaller amount of ra- 
dioactive yield and a rapid decay (timescales 
of days): coalescence of compact objects, ac- 
creting white dwarfs (0-Ne-Mg) and final He 
shell flash in AM CVn systems. 

Following Li & Paczynksi (1998), Kulka- 
rni (2005) considers the possibility of the de- 
bris of neutron star coalescence being heated 
by decaying neutrons. Amazingly (despite 
the 10-min decay time of free neutrons) such 
events (dubbed as "macronovae" ) are de- 
tectable in the nearby Universe over a pe- 
riod as long day, provided even a small 



■^This is a critical assumption and must be checked 
for every potential scenario under consideration. In 
a relativistic fireball most of the energy is transferred 
to matter. For novae, this assumption is violated 
(Shara, pers. comm.). 
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amount (~> 10~^ Mq) of free neutrons is re- 
leased in such explosions. 

Bildsten et al. (2007) consider a Helium 
nova (which arise in AM CVn systems). For 
these events (dubbed ".la" supernovae), not 
only radioactive Nickel but also radioactive 
Iron is expected. 

Intermediate mass stars present two possi- 
ble paths to sub-luminous supernovae. The 
0-Ne-Mg cores could either lead to a dis- 
ruption (bright SN but no remnant) or 
a sub-luminous explosion (Kitaura et al. 
2006). Separately, the issue of 0-Ne-Mg 
white dwarfs accreting matter from a com- 
panion continues to fascinate astronomers. 
The likely possibility is a neutron star but the 
outcome depends severely on the unknown ef- 
fects of rotation and magnetic fields. One 
possibility is an explosion with low Nickel 
yield (see Metzger et al. 2008 for a recent 
discussion and review of the literature). 

An entirely different class of explosive 
events are expected to arise in massive or 
large stars: birth of black holes (which can 
range from very silent events to GRBs and ev- 
erything in between), strong shocks in super- 
giants (van den Heuvel 2008) and common 
envelope mergers. Equations [5] and [6] pro- 
vide a guidance to the expected appearance 
of such objects. Fryer et al. 2007 develop 
detailed model for faint, fast supernovae due 
to Nickel "fallback" into the black hole. For 
the case of the birth of a black hole with 
no resulting radioactive yield (the newly syn- 
thesized material could be advected into the 
black hole) the star will slowly fade away 
on a timescale of min(rrf, r). Modern sur- 
veys are capable of finding such wimpy events 
(Kochanek et al. 2008). 

In the spirit of this open-ended discussion 
of new transients we also consider the case 
where the gas pressure could dominate over 
photon pressure. This is the regime of weak 
explosions. If so, P = 2/3{E/V) and Equa- 



tion [T] can be integrated to yield 
Lit) = ^ n • 

[t/Th + 1) V ThTd ) 

In this case the relevant timescale is the hy- 
drodynamic timescale. This regime is popu- 
lated by Luminous Blue Variables and Hyper- 
giants. Some of these stars are barely bound 
and suffer from bouts of unstable mass loss 
and photometric instabilities. 

As can be gathered from Figure [1] the pace 
of discoveries over the past two years gives 
great confidence to our expectation of filling 
in the phase space of explosions. 

3 Rationale & Motiva- 
tion: New Astronomy 

Four entirely different fields of astronomy are 
now being nurtured by physicists. For these 
fields all that matters is transients in the local 
Universe. 

Cosmic rays with energies exceeding 
10^° eV are strongly attenuated owing to 
the production of pions through interac- 
tion with the cosmic micro-wave back- 
ground (CMB) photons [the famous Greisen- 
Zatsepin-Kuzmin (GZK) effect]. Recently, 
the Pierre Auger Observatory (PAO 2007) 
has found evidence showing that such cosmic 
rays with energies above 6 x 10^^ eV are corre- 
lated with distribution of matter in the local 
75-Mpc sphere. 

Similarly very high energy (VHE) pho- 
tons (TeV and PeV) have a highly restricted 
horizon. The TeV photons interact with 
CMB photons and produces electron-positron 
pairs. A number of facilities are now rou- 
tinely detecting extra-galactic TeV photons 
from objects in the nearby Universe (Veritas, 
MAGIC, HESS, CANGAROO). 

Neutrino astronomy is another budding 
field with an expected vast increase in sensi- 
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tivity. The horizon here is primarily hmited 
by sensitivity of the telescopes (ICECUBE). 

Gravitational wave (GW) astronomy suf- 
fers from both poor localization (small in- 
terferometer baselines) and sensitivity. The 
horizon radius is 50 Mpc for enhanced LIGO 
(e-LIGO) and about 200 Mpc for advanced 
LIGO (a-LIGO) to observe neutron star co- 
alescence. There is wide spread agreement 
that the greatest gains will require electro- 
magnetic localization of the event. 

Progress in these areas, especially GW as- 
tronomy, require arc-second localization. Un- 
fortunately, all sorts of foreground and back- 
ground transients will be found within the 
several to tens of sq deg of expected localiza- 
tions. Studying each of these transients will 
result in significant "opportunity cost" . This 
issue is addressed in ^ 

4 Foreground Fog & 
Background Haze 

Ongoing projects of modest scope offer a 
glimpse of the pitfalls on the road to achiev- 
ing the grand goals summarized in [|2] and ^ 
Nightly monitoring of M31 for novae (several 
groups) and our Palomar 60-inch program 
of nearby galaxies already show high vari- 
ance in the MMRD relation (Figure [2]). The 
large scatter of the new novae suggests that 
in addition to the mass of the white dwarf, 
other physical parameters play a role (such 
as accretion rate, white dwarf luminosity, e.g. 
Shara 1981). 

A nightly targeted search of nearby rich 
clusters (Virgo, Coma and Fornax) using the 
CFHT (dubbed "COVET") and the 100-inch 
du Pont (Rau et al. 2008a) telescopes uncover 
the extensive fore-ground fog (asteroids, M 
dwarf flares, dwarf novae) and the back- 
ground haze (distant, unrelated SN). The pie- 
chart in Figure [3] dramatically illustrates that 



new discoveries require efficient elimination 
of fore- and hack- ground events. 
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Fi gUre 2) A plot of the peak absolute magnitudes ver- 
sus decay timescale of novae discovered by the Palomar P60- 
FasTING project. The shaded grey region represents the Max- 
imum Magnitude Rate of Decline relationship bounded by 
±3(7 (Delia Valle and Livio 1995). The data that defined this 
MMRD are shown by green circles. Squares indicate novae dis- 
covered by P60-FasTING in 2007-2008 (unlabelled if in M31). 
[Preliminary results from Kasliwal et al. 2009b, in prep.] 




Figure 3: 28 covet transients were discovered dur- 
ing our pilot run in 2008A (7 hours) — two novae and the re- 
mainder background supernovae and AGN. Transients with no 
point source or galaxy host to a limiting magnitude of r>24 
are classified as hostless. Of the 2800 candidates, our pipeline 
automatically rejected 99% as solar-system or galactic. [Pre- 
liminary version from Kasliwal et al. 2009c, in prep.] 

5 The Era of Synoptic 
Imaging Facilities 

There is widespread agreement that we are 
now on the threshold of the era of synop- 
tic and wide field imaging at optical wave- 
lengths. This is best illustrated by the profu- 
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sion of operational (Palomar Transient Fac- 
tory, PanSTARRS), imminent (SkyMapper, 
VST, ODI) and future facilities (LSST) in op- 
tical and SASIR in infrared. 

In Table [1] and Figure H] we present current 
best estimates for the rates of various events 
and the "grasp" of different surveys. 




-20 -18 -16 -14 -12 

Transient Peak M (mag) 



Fi gUre 4; Volume probed by various surveys (in a single 
specified cadence period) as a function of transient absolute 
magnitude. Red plusses represent the minimum survey vol- 
ume needed to detect a single transient event (the uncertainty 
in the y-axis is due to uncertainty in rates). PTF-5d (blue- 
solid) is more sensitive than TSS (dotted), Skymapper (dot- 
dashed), Supernova Legacy Survey (SNLS, dashed) and SDSS- 
SN (double-dot dashed) and competitive with PanSTARRS- 1 
(PSl-MD, long dashed). Lines for each survey represent one 
transient event in specified cadence period. PTF-ld (green 
solid line) represents a targeted 800 sq dcg survey probing lu- 
minosity concentrations in the local Universe, with a factor of 
three larger effective survey volume than a blind survey with 
same solid angle. PTF will discover hundreds of supernovae 
and possibly several rare events such as "O.Ia", Luminous Red 
Novae (LRNe) and Macronovae (MNe) per year. The LSST 
(Wide Fast Survey) will discover hundreds of these rare events. 
[Adapted from figure by Bildstcn ct al. 2009.] 

The reader should be cautioned that many 
of these rates are very rough. Indeed, the 
principal goal of the Palomar Transient Fac- 
tory is to accurately establish the rates of the 
fore- and back-ground events. It is clear from 
Figure m that the impressive grasp of LSST is 
essential to the discovery of rare events. 



6 Localizing LIGO events 
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Fi gUre 5 1 Number of false positives in a single snapshot. 
For e-LIGO (blue circles), we use the median localization of 41 
sq deg and follow-up depth of r<21 is assumed. For a-LIGO 
(red stars), we use the median localization of 12 sq deg and 
follow-up depth of r<24. Filled symbols denote false positives 
in the entire error circle and empty symbols show false positives 
that are spatially coincident with nearby galaxies. Dwarf novae 
and M-dwarf (dM) flares constitute the foreground fog and the 
"error bars" on numbers represents the dependence on galactic 
latitude. Supernovae (Ia,IIp) constitute background haze. 



act localization on sky (assuming a neutron- 
star neutron-star merger waveform and triple 
coincidence data from LIGO-Hanford, LIGO- 
Louisiana and Virgo). The localizations 
range between 3-700 sq deg for e-LIGO and 
0.2-300 sq deg for a-LIGO (range quoted be- 
tween 10th and 90th percentile). The Uni- 
verse is very dynamic and the number of false 
positives in a single snapshot is several tens 
for a median localization (see Figure E]). 

Fortunately, the sensitivity-limited < 
200 Mpc horizon of GW astronomy is a bless- 
ing in disguise. The opportunity cost (§!]) can 
be substantially reduced by restricting follow- 
up to those transients that are spatially coin- 
cident with galaxies within 200 Mpc. Limit- 
ing the search to the area covered by galax- 
ies within a LIGO localization reduces a sq 
deg problem to a sq arcmin problem — a re- 
duction in false positives by three orders of 
magnitude! 



We simulated a hundred GW events (Kasli- Given the total galaxy light in the local- 
wal et al. 2009a, in prep) and compute the ex- ization, we also find that the number of false 
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Table 1: Properties and Rates for Optical Transients" 



Class 


M„ 




T 


Universal Kate (UKj 


it Kate 


T COT* ~D ^4-^ 

Lbbi Kate 




[mag] 


[days] 






[yr-i] 


Luminous red novae 


-9.. - 


13 


20. .60 


(1..10) X 10-1^ yr-i L-V 


0.5. .8 


80.. 3400 


Fallback SNe 


-4.. - 


21 


0.5. .2 


< 5 X 10-6 Mpc-3 yr-'i 


<3 


<800 


Macronovae 


-13.. - 


15 


0.3. .3 


lO-^' -^Mpc-^ yj.-l 


0.3. .3 


120. .1200 


SNe .la 


-15.. - 


17 


2. .5 


(0.6. .2) X lO-^Mpc-^ yr-i 


4.. 25 


1400.. 8000 


SNe la 


-17.. - 


19.5 


30.. 70 


3 X lO-^Mpc-^ yr-i 


700 


200000^* 


SNe II 


-15.. - 


20 


20. .300 


(3. .8) X lO^^Mpc-^ yr-i 


300 


lOOOOO'' 



"Table from Rau et al. 2008b; see references therein. ''Time to decay by 2 magnitudes from peak. 
■^Universal rate at z< 0.12. ''From M. Wood-Vassey, pers. comm. 



Table 2: Galaxy Characteristics in LIGO Localizations 







E-LIGO 






A-LIGO 






10% 


50% 


90% 


10% 


50% 


90% 


GW Localization (sq deg) 


3 


41 


713 


0.2 


12 


319 


Galaxy Area (sq arcmin) 


4.4 


26 


487 


0.15 


20.1 


185 


Galaxy Number 


1 


31 


231 


1 


76 


676 


Galaxy Luminosity (Log) 


10.3 


11.2 


12.1 


10.9 


12.0 


13.0 



positives due to unrelated supernovae or no- 
vae within the galaxy is negligible. To be 
sensitive to transients as faint at peak as 
M<— 13 (fainter than the faintest observed 
short hard gamma ray burst optical after- 
glow), e-LIGO needs at-least a 1-m class tele- 
scope for follow-up (m<21, 50Mpc) and a- 
LIGO an 8-m class (m< 24, 200 Mpc). Given 
the large numbers of galaxies within the local- 
ization (Table [2]), a large field of view camera 
(> 5 sq deg) will help maximize depth and 
cadence as compared to individual pointings. 
Thus, PTF is well-positioned to follow-up e- 
LIGO events and LSST to follow-up a-LIGO 
events. 

The feasibility hinges on a complete galaxy 
catalog. We compiled all available distances 
to nearby galaxies and find that this catalog 
is only 70% complete at e-LIGO distance 
and 55% complete at a-LIGO distance. It 
is important that astronomers embark upon 
completing the census of galaxies within 
200 Mpc. 
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